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1 introduction 

In this note I announce and introduce the program LCFIT developed for fitting harmonic 
functions to a data set, particularly to time-series data. LCfit stands for Linear Combi- 
nation fitting. 

The program has linear and non-linear fitting modes and uses an imple mentation 



of the Levenberg-Marquardt (LM) least squares fitting algorithm written by lLourakis 



(120041 ) for optimizing the parameters. The user gives a data set, base frequencies and 
linear-combination terms as input. In linear mode, LCfit fits the phases and amplitudes 
of the linear-combination terms to the data. In non-linear mode, the base frequencies are 
also fitted using the user input as initial values for these parameters. The result of the fit 
is given in different forms. The formal uncertainties of the fitted parameters, calculated 
from the covariance matrix of the results are also given in the output. The program 
is written in C++. The source code is freely distributed under the terms of the GNU 
General Public Licens^] and is available from the authorj^' 



2 usage 

Brief instructions on the usage is written on the standard output if the program is run 
without sufficient number of command-line arguments: 



J http : //www . gnu . org/ copylef t/gpl . html 

2 Web page: http://konkoly.hu/staff/sodor/lcfit.html 

3 e-mail: sodor *at* konkoly *dot* hu 
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Usage: lcfit [options] data fitdef [pw_output] 



where : 
options: 



data 

fitdef 

pw_output 



-g 



-V 



-c {file} 



gnuplot output 
verbose 

mufran . coe output to file {file} (if used, -c option must be the last one) 

file containing the light curve (time series) 

file describing the frequencies and lincombs to fit. 

file to write prewhitened output. Hill be overwritten! 



The data file is a simple list of whitespace separated time value pairs. Empty lines 
and lines beginning with '#' are ignored. The fitdef file format is introduced in the next 
section. If the pw_output command-line argument is given, then a time series prewhitened 
with the fitted solution is written into that file. 



The harmonic functions to fit, as well as some parameters of the fitting algorithm are 
given in the fitdef input file, which has the following syntax: 

[f ixf req] 

sid = 1.0024 

[varf req] 

fl = 2.1453 
f2 = 0.0149 

[lincomb] 

sid 
f2 

f 1 

2*f 1 

fl + f2 
3*fl - 2*f2 

[LMpar] 

maxiter = 5000 

#epsilon2 = 3e-7 

#mu=le-8 

tO = 54540.5 



The fitdef file is divided into four sections by the headings: [f ixf req] , [varf req] , 
[lincomb] , and [LMpar] . The section headings are case sensitive and their order is 
fixed. However, not each of them is mandatory. One of the frequency-definition sections 
( [f ixf req] or [varf req]) and [LMpar] might be omitted, but one of the frequency 
definition sections and [lincomb] are mandatory. Empty lines and lines beginning with 



2.1 The fitdef file 
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'#' mark are ignored, as well as whitespaces around the +, -, * and = marks and at the 
beginning and end of lines. 

The base frequencies not to be fitted are given in the [fixfreq] section in 
identifier = value format, one on each line. The frequency unit is the inverse of 
the unit of time used in the input data time-series file. The frequencies given in the 
[varf req] section will be fitted by LCFIT. It is mandatory to give initial frequency val- 
ues for the fitted frequencies. If at least one frequency is given in the [varf req] section, 
a non-linear fit will be calculated. 

The actual Fourier terms to be fitted are given in the [lincomb] section. Note that 
harmonic terms of the base frequencies will be fitted only when the given base frequency 
is listed here. 

The last, optional section, [LMpar] lists further parameters for the program: tO is 
the epoch for which the output should correspond (its default value is the integer part 
of the time of the first data point of the input data file). The maximum number of 
iterations of the LM algorithm is given by maxiter (its default value is 1500, which might 
be insufficient for a complex fit). The step-size for numeric difference calculations, which 
is set to its default by LCFIT according to the length of the data set in order to fit 
the frequencies appropriately, is given by delta. The modification of this parameter by 
hand is usually not necessary. For completeness, we mention mu, epsilonl, epsilon2, 
and epsilon3, which are further LM p arameters (for details we refer the reader to the 



documentation of the LM algorithm in iLourakid 12004 ) . Their default values are almost 



always suitable, thus, it is seldom necessary to override them. 
2.2 Output 

The result of the fit is written to the standard output in the following format: 

Final parameters: 



Results for harmonic curve fit: 

LM stopped after 70 iterations and 95 function evaluations. 
Exit reason: convergence (2 - stopped by small Dp) 

Minimization info: 

Final r.m.s. = 0.00984932 Reduced Final r.m.s. = 0.00986411 



sid = 1.0024000000; # (fix freq.) 
fl = 2.1452983111; Dfl = 1.0000000000 

f2 = 0.0148972133; Df2 = 0.0000020287 

0.00011 0.00031 55490.00000 # AO DA0 epoch 



#N Al 



001 
002 
003 
004 
005 
006 



0.001546 
-0.003992 
-0.473363 
-0.247322 
-0.007902 
-0.000214 



DAI 

0.000002 
0.000394 
. 000497 
0.000640 
0.913895 
0.540259 



A2 



DA2 



1.001022 0.000442 
0.000678 0.000040 
-0.123707 0.000090 
0.056327 0.000077 
-0.003489 0.403028 
-0.000333 0.841428 



A DA f i 

1.001023 0.000441 1. 

0.004049 0.000376 2. 

0.489261 0.000570 3. 

0.253655 0.000599 2. 

0.008638 0.000546 3. 

0.000396 0.000654 4. 



Dfi 

569252 0.000885 
973280 0.103134 
397212 0.001550 
917664 0.003001 
557396 0.254601 
141593 1.124066 



1.0024000000 
0.0148972133 
2.1452983111 
4.2905966221 
2.1601955244 



lincomb 
(+l*sid) 
(+I*f2) 
(+l*f 1) 
(+2*f 1) 
(+l*f 1+I*f2) 



6.4061005065 (+3*f l+-2*f 2) 



In the first part, a summary can be read about the outcome of the fit. 
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The Minimization info section is self-explanatory. 

The next section gives the results of the frequency fit. The example shows the output 
of a fit with both fixed (sid) and fitted (f 1, f2) freqencies, according to the example 
f itdef file. The formal uncertainties are also given for the fitted frequencies. 

After the frequency section follows a single line with the fitted zero-point of the data 
(AO), its uncertainty (DAO) and the epoch or to to which the Fourier parameters listed in 
the last section correspond. 

The last section lists the Fourier parameters of the fitted linear-combination terms. 
The first column is a serial number of the term, followed by Al and A2, the amplitudes 
in cosine- and sine-term (or real and imaginary part) formalism, while A and fi are 
amplitudes and phases in the sine-amplitude and -phase formalism. The amplitudes are 
given in the same unit as used in the input file, while the phases are in radians. Dx denotes 
the uncertainty of x. The last two columns give the frequency of the term (f ) and the 
linear-combination expression (lincomb). 

In case the -g (GNUPLOT output) option is selected on the command line, then the 
results written to the standard output are formatted in a different way so that the output 
can directly be used as input for GNUPLOT as in the following example: 



$ lcfit -g lightcurve_f ile fitdef_file I gnuplot -persist 



As a result, gnuplot shows the data and the fitted curve for inspection. 

An alternative output format can also be chosen with the -c {file} command line 
option. In this case, the results are not only printed on the standard output, but the 
Fourier par ameters are al so written into {file} in the Fourier coefficient file format of 
MuFrAn (lKollathl[l99nh . 

In verbose mode, further information are printed during the run of LCFIT both on 
the standard error and on the standard output channels. If the pw_output command-line 
argument is given, then a time series prewhitened with the fitted solution is written into 
that file. 
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